Popular Summary 


What would have happened to the ozone layer if chlorofluorocarbons (CFCs) had 

not been regulated? 


In 1974, Drs. Mario Molina and F. Sherwood Rowland published a landmark paper that 
suggested that ozone depletion could be caused by human produced compounds known 
as chlorofluorocarbons (CFCs). Ozone is the Earth’s natural sunscreen, so depletion of 
the ozone layer would lead to an increase in skin cancer causing ultraviolet radiation. 
Since 1974, the scientific connection between ozone losses and CFCs and other ozone 
depleting substances (ODSs) has been firmly established with laboratory measurements, 
atmospheric observations, and modeling research. The severe depletion of ozone over 
Antarctica, the so-called ozone hole, is a result of the release of chlorine from these CFC 
compounds. 

This science research led to the implementation of the Montreal Protocol in 1987 by the 
countries of the world. The Protocol and further amendments largely stopped the 
production of ODSs by 1992, and the levels of these gases are now declining in both the 
lower and upper atmosphere. It is expected that the ozone layer should be recovered in 
the mid-latitudes by 2050 and in the polar latitudes by 2070. In 1995, Drs. Mario Molina 
and F. Sherwood Rowland were awarded the Nobel Prize in Chemistry for their work on 
ozone depletion. 

What would have happened to our atmosphere if CFCs had never been regulated? We 
have used a fully coupled radiation-chemical-dynamical model to simulate a future world 
where ODS production grew at an annual rate of 3%. By 2065, 67% of the overhead 
ozone is destroyed in comparison to 1980. Large ozone depletions in the polar region 
become year-round rather than just seasonal as is currently observed in the Antarctic 
ozone hole. Ozone levels in the tropical lower stratosphere remain constant until about 
2053 and then collapse to near zero by 2058 as a result of “polar ozone hole” chemical 
processes developing in the tropics. In response to ozone changes, ultraviolet (UV) 
radiation increases, more than doubling the “sun-burning” radiation in the northern 
summer mid-latitudes by 2060. 
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Abstract. Ozone depletion by chlorofluorocarbons (CFCs) was first proposed by Molina and Row- 
land in their 1974 Nature paper. Since that time, the scientific connection between ozone losses 
and CFCs and other ozone depleting substances (ODSs) has been firmly established with labora- 
tory measurements, atmospheric observations, and modeling research. This science research led to 
5 the implementation of international agreements that largely stopped the production of ODSs. In 
this study we use a fully-coupled radiation-chemical-dynamical model to simulate a future world 
where ODSs were never regulated and ODS production grew at an annual rate of 3%. In this “world 
avoided” simulation, 17% of the globally-average column ozone is destroyed by 2020, and 67% 
is destroyed by 2065 in comparison to 1980. Large ozone depletions in the polar region become 
1 0 year-round rather than just seasonal as is currently observed in the Antarctic ozone hole. Very large 
temperature decreases are observed in response to circulation changes and decreased shortwave ra- 
diation absorption by ozone. Ozone levels in the tropical lower stratosphere remain constant until 
about 2053 and then collapse to near zero by 2058 as a result of heterogeneous chemical processes 
(as currently observed in the Antarctic ozone hole). The tropical cooling that triggers the ozone 
1 5 collapse is caused by an increase of the tropical upwelling. In response to ozone changes, ultraviolet 
radiation increases, more than doubling the erythemal radiation in the northern summer midlatitudes 
by 2060. 
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1 Introduction 


Molina and Rowland (1974) were the first to propose that ozone could be depleted by the release 
20 of chlorofiuorocarbons (CFCs) to the atmosphere. The chemical breakdown of CFCs and other 
ozone depleting substances (ODSs) in the stratosphere releases chlorine (Cl) and bromine (Br) atoms 
that destroy ozone molecules in catalytic cycles. Ozone depletion would result in an increase of 
biologically harmful solar ultraviolet (UV) radiation. Following this discovery, public reactions to 
potential ozone depletion led to a ban on the use of CFCs as aerosol propellants by the late 1970s, 
25 slowing CFC growth rates. However, by 1980 CFC production had again begun to rise as their use 
as blowing agents and refrigerants grew. Continued scientific research led to the writing of a series 
of reports that culminated in the first World Meteorological Organization (WMO)/United Nations 
Environment Program scientific assessment of ozone depletion (WMO, 1985). This assessment 
provided the scientific consensus that CFCs and halons posed a serious threat to the ozone layer. In 
30 response, the landmark Montreal Protocol agreement was negotiated in 1 987 (Sarma and Bankobeza, 
2000). This agreement regulated the production of chlorofluorocarbons and other ODSs. Since 1987, 
193 nations have signed the Montreal Protocol . Five amendments to the Montreal Protocol have now 
led to the cessation of the major portions of ODS production around the world. The cumulative levels 
of chlorine and bromine from the ODSs are now decreasing in both the troposphere (Montzka et al., 
35 1996) and the stratosphere (Anderson et al., 2000; Froidevaux, et al., 2006). The work of Molina and 

Rowland (1974) led to a major reorganization and redirection of chemical production to substitute 
compounds, to changes in technologies that traditionally used CFCs, and to the first international 
agreement on regulation of chemical pollutants. 

The regulation of ODSs was based upon the ozone assessments that presented the consensus 
40 of the science community. The regulation presupposed that a lack of action would lead to severe 
ozone depletion with consequent severe increases of solar UV radiation levels at the Earth’s surface. 
Because of the successful regulation of ODSs, ozone science has now entered into the accountability 
phase. There are two relevant questions in this phase. First, are ODS levels decreasing and ozone 
increasing as expected because of the Montreal Protocol regulations? Second, what would have 
45 happened to the atmosphere if no actions had been taken? That is to say, what kind of world was 
avoided? It is this latter question that is the focus of the present study. 

Two studies have modeled the impact of high levels of CFCs on ozone. First, Prather et al. (1996) 
modeled the ozone response to continued growth of the ODSs without the Montreal Protocol and 
its amendments. They hypothesized that ODSs continued to grow at “business-as-usual” growth 
50 rates from 1974 into the early 21st century. Such ODS growth could have resulted from a lack of 
understanding of ozone depletion, a failure to regulate ODS production, and further development of 
technologies that used CFCs. In their study, ODSs were allowed to grow by approximately 5-7% 
per year, resulting in a total stratospheric chlorine (Cl y ) level of about 9ppbv by 2002. They used a 
two-dimensional (2-D, or latitude-altitude) model with projected CFC and halon levels to calculate a 


2 



55 globally-averaged total ozone depletion of 10% by the year 1999. Second, Morgenstern et al. (2008) 
used a three-dimensional (3-D, or longitude-latitude-altitude) coupled chemistry-climate model to 
simulate the impact of 9ppbv on the stratosphere in a timeslice or fixed time mode (equivalent to 
levels of Cl y that they estimated would have been reached in about 2030). In their simulation they 
found peak ozone depletions (> 35%) in the upper stratosphere, with additional ozone depletion in 
60 both polar regions (> 20%). 

In addition to the ozone studies, Velders et al. (2007) have recently investigated the impact of 
unrestrained growth of CFCs on climate. They pointed out that because CFCs and other ODSs are 
also greenhouse gases, the ODSs would have added an additional 0.8-1 .6 W m~* 2 of radiative forcing 
by 2010 if the ODSs had continued to grow at 3-7% per year after 1974. Hence, the Implementation 
65 of the Montreal Protocol has limited both ozone depletion and climate change. 

Prather et al. (1996) was limited by using a fixed transport 2-D model (Jackman, et al., 1996) 
to produce their estimates of ozone loss. Since this version of the 2-D model did not include the 
feedback among radiative, dynamical, and chemical processes, the model simulations were termi- 
nated when total ozone depletion reached 10%. Such significant computed total ozone depletion 
70 was assumed to cause large stratospheric changes, both radiatively and dynamically. Models of the 
stratosphere have greatly improved over the last decade. Using 3-D instead of 2-D models bet- 
ter account for longitudinal variations and more realistically represent the polar latitude dynamics. 
State-of-the-art models now have complex feedback among the chemical, radiative, and dynamical 
processes. 

75 Here we employ the Goddard Earth Observing System (GEOS) chemistry-climate model (CCM) 
to explore the ozone distributions that might have been without ODS regulations (a “world avoided”), 
updating and extending the results of Prather et al. (1996) up to 2065. Because ozone is a radiatively 
active gas in both the ultraviolet and the infrared, the destruction of ozone alters the temperature 
and the wind distributions, which then affects the transport of ozone and other gases. Further, as 
80 temperatures change, ozone loss rates change, producing a feedback that further alters the ozone 
distribution. Hence, a proper simulation of the “world avoided” requires a CCM with interactive 
radiation, chemistry, and dynamics. 

This “world avoided” estimate of ozone depletion is performed for a few basic reasons. First, 
scientists predicted massive ozone losses, actions were taken, and these losses have not occurred. 
85 Hence, do the state-of-the-art models actually predict the large losses that were hypothesized in the 
1980s? Second, such model estimates give an indication of the unforeseen impacts of large ozone 
losses on the chemistry and climate of the stratosphere. Such extreme simulations are useful for 
evaluating theoretical expectations of dynamics and chemistry (e.g., how does the tropopause change 
if the ozone layer is destroyed?). Third, a transient simulation provides an ensemble of total chlorine 
90 and bromine values and their associated ozone losses, from small chlorine values in the 1960s to 
high chlorine values in the 2060s. Because we cannot actually predict how large CFC and other 
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ODS concentrations might have become, this simulation provides a large range of CFC levels from 
which to choose. Fourth, this simulation provides a diagnostic test of the model, pushing it into an 
extreme chemical-radiative-dynamical state that could uncover subtle non-physical model problems. 

95 Finally, and most importantly, this type of estimate provides a quantitative baseline for assessing 
the impact of international agreements on ozone, UV radiation impacts, and climate change. The 
acknowledged success of the Montreal Protocol in the present must be measured against what might 
have been without that agreement. 

The scope of this study is limited to describing the very large ozone losses and the subsequent 
1 00 rise of ultraviolet (UV) levels. We only include highlights of the radiative, chemical, and dynamical 
effects in the stratosphere. There are eight sections including the introduction. Section 2 describes 
the experiment with subsections on the model and the details of the simulations. Section 3 describes 
the evolution of halogens in the model experiments. Section 4 describes the evolution and distri- 
butions of ozone. Sections 5 and 6 discuss some important aspects of the chemistry and dynamics, 
105 respectively. Section 7 provides estimates of surface ultraviolet levels. The final section provides a 
summary and discussion of some important aspects of this “world avoided” scenario. 

2 Experiment 

The model simulations use the GEOSCCM, which is described in detail by Pavvson et al. (2008). 

2.1 Model description 

1 1 0 Briefly, the model has a horizontal resolution of 2° latitude by 2.5° longitude with 55 vertical levels 
up to 0.01 hPa (80km). The dynamical time step is 7.5 minutes. The model uses a flux-form semi- 
Lagrangian dynamical core (Lin, 2004) and a mountain-forced gravity- wave drag scheme with a 
wave forcing for non-zero phase speeds (Garcia and Boville, 1994; Kiehl, et al., 1998). The sub-grid 
moisture physics and radiation are adapted from Kiehl, et al. (1998), as described by Bosilovich 
115 et al. (2005). The photochemistry code is based on the family approach, as described by Douglass 
and Kawa (1999) and uses the chemical kinetics from Jet Propulsion Laboratory Evaluation 14 
(Sander et al., 2002). Tropospheric ozone is relaxed to a climatology (Logan, 1999) with a 5-day 
time scale. Processes involving polar stratospheric clouds (PSCs)use the parameterization described 
by Considine et al. (2000). A total of 35 trace gases are transported in the model. Gases with 
120 surface sources are specified as mixing ratios in the lowest model layer (e.g., “greenhouse gases” 
such as CO2 and “ozone-depleting gases” such as CFC-1 1). The coupling between the chemical 
and physical state of the atmosphere occurs directly through the GCM’s radiation code. The model- 
simulated H 2 0, C0 2 , 0 3 , CH 4 , N 2 0, CFC-1 1, CFC- 12, and HCFC-22 are used in the radiation 
computations. The same Alb Intergovernmental Panel on Climate Change (IPCC) greenhouse gas 
1 25 scenario is used for C0 2 , CH4 , and N 2 0 in all of our simulations (Nakicenovic and Swart, 2000). 
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GEOSCCM has been compared to observations and other models and does very well in repro- 
ducing past climate and key transport processes (Eyring, et al., 2006; Stolarski et al., 2006; Oman 
et al., 2008; Pawson et al., 2008). Pawson et al. (2008) showed good agreement between the ozone 
simulations and observations, but noted a slightly high total ozone bias and a too cold and long-lived 
1 30 Antarctic vortex. Predictions from the GEOSCCM have been compared to other model predictions 
(Eyring, et al., 2007), showing good general agreement with other CCMs. 

2.2 Model simulations 

There are four simulations used in this study: a “ reference past”, a “reference future ”, a “ fixed chlo- 
rine” (with chlorine levels fixed at the 1960 level), and the “ WORLD AVOIDED”. Table 1 summa- 
1 35 rizes the four simulations. Figures in this study that compare the simulations use color to distinguish 
between them, using blue for the reference past , red for the reference future, green for the fixed 
chlorine , and black for the WORLD AVOIDED. The first three simulations have been described in 
Oman et al. (2008) and Pawson et al. (2008). The reference past simulations have also been shown 
by Eyring, et al. (2006), while the reference future also appears in Eyring, et al. (2007). 

140 The reference past simulation is an attempt to simulate the past observations of the atmosphere. 
Hence this simulation is forced by the observed ODSs at the surface (Montzka et al., 2003), and 
uses the observed sea surface temperatures (SSTs) and ice distribution (Rayner et al., 2003). The 
reference future is an attempt to simulate the future using our best guesses for future ODSs (Montzka 
et al., 2003) and SSTs. There are two reference future simulations shown herein. All of the reference 
1 45 future plots versus time use the same simulation shown by Eyring, et al. (2007) that was forced with 
the HadGEMl SSTs. All of the zonal mean difference plots use a reference future simulation that 
is forced with the National Center for Atmospheric Research (NCAR) Community Climate System 
Model version 3 (CCSM3) SSTs (Collins et al., 2006). This second reference future simulation (B 
in Table 1) is used because the SSTs are consistent with the WORLD AVOIDED simulation in the 
150 period from 2050-2065. The differences between these two reference future simulations are quite 
small in the stratosphere. 

The WORLD AVOIDED simulation extends from 1974 through 2065. This simulation is driven 
by one of the mixing ratio scenarios established in Velders et al. (2007), where ODS estimates are 
based upon a scenario (their MR74) wherein production of CFCs and halons grows annually at 3% 
155 beginning in 1974. In this scenario, there was no early warning of the danger of CFCs as actually 
occurred because of Molina and Rowland (1974). The 3% per year production growth is lower than 
the observed growth of 12-17% in CFC production in the period up to 1974. The surface total 
chlorine reaches a level of 9 ppbv in 2012. In contrast, the CFC growth used in Prather et al. (1996) 
was approximately 7% and Cl y reached a level of 9 ppbv in about 2002. Morgenstern et al. (2008) 
1 60 compared two fixed Cl y levels of 3.5 ppbv and 9 ppbv in their timeslice simulation. In our WORLD 
AVOIDED simulation, the upper stratospheric Ci y reaches 9 ppbv in about 2019. 
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Reliable estimates of surface temperature changes can only be made using coupled atmosphere- 
ocean general circulation models (AOGCMs) (Houghton, et al., 1990). In our CCM we attempt to 
overcome this problem by using SSTs that are consistent with increasing greenhouse gases. The 
165 SSTs and sea ice used in the WORLD AVOIDED run are derived from two integrations of the 
NCAR Community Climate System Model, from the version 2 (CCSM2) for 1974-2049 and from 
the CCSM3 for 2050-2065. The latter is from an IPCC AR4 integration, but both are consistent with 
the IPCC Alb scenario growth of greenhouse gases. A limitation of this prescribed SST approach 
is that it is inconsistent with the growth of the ODSs in our WORLD AVOIDED simulation. Since 
1 70 ODSs are also greenhouse gases, they introduce a direct radiative forcing in the atmosphere (less so 
in the stratosphere since CFCs are destroyed there). The prescribed SST used herein acts as a thermal 
sink, constraining the troposphere to evolve in parallel with Alb greenhouse gas growth rather than 
with the additional radiative forcing from the higher levels of CFCs. Assuming a climate sensitivity 
parameter of 0.5 K (W m” 2 )"' 1 (Ramanathan et al., 1985), the additional CFCs would contribute a 
175 surface warming of very roughly 0.25 K by 2010 following the Velders et al. (2007) estimates of 
the CFC induced radiative forcing. Hence, the WORLD AVOIDED troposphere and oceans are too 
cold compared with what would develop for the same scenario using an AOGCM. Furthermore, the 
extreme UV levels simulated here would also impact tropospheric chemistry, tropospheric ozone 
radiative forcing, and biogeochemical processes. The GEOSCCM model is a stratosphere-only 
180 chemistry model that uses a climatological distribution of ozone in the troposphere with a 5-day 
relaxation time scale. While tropospheric reductions of ozone are evident in the model’s ozone field, 
these reductions are solely a result of the advection of low ozone from the stratosphere. 

We also simulated the WORLD AVOIDED using a 2-D model with the Velders et al. (2007) MR74 
scenario. This 2-D model (Rosenfield et al., 2002) is a newer version of the 2-D model used by 
185 Prather et al. (1996). This newer model has coupled chemistry, radiation, and dynamics. The dy- 
namics of the newer 2-D model has a parameterization of eddy effects that interact with the mean 
circulation. While the dynamics, chemistry, and physics of the GEOSCCM and the 2-D model are 
very distinct, the newer 2-D model results are in good quantitative agreement with the results herein 
from the GEOSCCM WORLD AVOIDED simulation. 

1 90 3 Atmospheric halogens 

In this WORLD AVOIDED simulation, stratospheric ODS levels increase rapidly after 2000. Figure 1 
shows the equivalent effective chlorine (EEC1) from Velders et al. (2007) (grey line) that is used as 
the surface mixing ratio boundary condition in the WORLD AVOIDED simulation. This EEC! is 
estimated by totaling all of the inorganic chlorine and bromine separately and then adding them with 
195 the inorganic bromine multiplied by a factor of 60 to account for the greater efficiency of bromine 
for ozone destruction. 
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The effective equivalent stratospheric chlorine (EESC) is calculated in Fig.l from the model out- 
put by summing the globally-averaged Cl y and Br y at 4.5 hPa with a Br y scaling factor of 60 (thick 
black line). We also theoretically estimate EESC following Newman et al. (2007), using the surface 
200 mixing ratio estimates and fractional release rates for each species, and an age spectrum with a 6-year 
mean age. This theoretical EESC (not shown here) overlaps the model EESC, giving us good confi- 
dence in these model estimates. In addition, Froidevaux, et al. (2006) estimated total chlorine from 
Microwave Limb Sounder EC1 observations in the upper stratosphere as approximately 3.6 ppbv in 
2006 with a slow decrease of about 0.8% per year. The reference future simulation has a peak in 
205 total chlorine of 3.4 ppbv with a decrease that matches the observed 0.8% per year, again giving us 
good confidence in our ability to simulate chlorine levels. The EESC theoretical estimates (again 
following Newman et al., 2007) are shown for the Montreal Protocol, the London Amendments, the 
Copenhagen Amendments, observations, and scenario Ab (Montzka et al., 2003). The current Al 
scenario (Daniel et al., 2007) is not shown here, but nearly overlaps the Ab scenario in Fig.l . 

210 The natural level of EESC is estimated to be approximately 1.2 ppbv (noted at the left in Fig. 1 . 
In the Ab scenario, the EESC reached a peak level of about 4.3 ppbv in approximately 2002. The 
WORLD AVOIDED scenario has EESC increasing to 1 1.5 ppbv by 2020 and 29.6 ppbv by 2050. In 
contrast, scenario Ab has EESC falling to 3.8 ppbv in 2020 and 2.8 ppbv by 2050. 

4 Ozone 

215 Global annually-averaged total ozone falls precipitously from about 3 15 DU in 1974 to about 1 10 DU 
in 2065 in the WORLD AVOIDED simulation. Approximately 27 DU of this annually-averaged total 
ozone is in the troposphere over the entire period. Figure 2 displays the total ozone levels from the 
four simulations. A value of less than 220 DU is nominally used to estimate the location and areal 
extent of the Antarctic ozone hole. The global annual average passes 220 DU shortly before 2040. 
220 In addition to the GEOSCCM WORLD AVOIDED simulation, we have also performed a WORLD 
AVOIDED simulation using the same ODS scenario with the coupled chemistry 2-D model (Rosen- 
held et al., 2002). The coupled 2-D model shows less ozone prior to 1980 and more global ozone in 
the 21st century with an offset of approximately 20 DU in 2040. Nevertheless, in spite of the large 
differences in these model implementations (e.g., 2-D vs. 3-D, chemical schemes, dynamical pa- 
225 rameterizations, etc.), the two models show good agreement over the simulation period and provide 
confidence in the magnitudes of the ozone losses. 

The sensitivity of total ozone to chlorine is roughly linear over the time period simulated. The 
grey-shaded inset to Fig. 2 shows the total ozone plotted against the model EESC. The ozone level 
is linearly anti -correlated with the EESC level. While the temporal tendency in Fig. 2 shows an 
230 accelerating decline (i.e., the change between 2040 and 2060 is greater than the change between 
2020 and 2040), this is a result of the accelerating increase of ODSs, rather than a greater sensitivity 
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to ODSs. 


Substantial reduction of ozone occurs at all latitudes in the WORLD AVOIDED simulation. Fig- 
ure 3 shows annually-averaged total ozone for a series of years from all four simulations. The 
235 total ozone in 1970 is represented by the fixed chlorine and reference past simulations (both have 
comparable EESC levels in 1970, see Fig. 1). These simulations are in reasonable agreement with 
observations (see Pawson et al., 2008) with a global total ozone average of about 310DU. The ref- 
erence future and WORLD AVOIDED simulations show substantial ozone losses by 2000 as ODSs 
increase. This loss trend continues as ODSs increase, such that global total ozone has fallen to a 
240 value of about 120 DU by 2060 in the WORLD AVOIDED simulation. The midlatitude maximum 
in both hemispheres has virtually disappeared by 2065, and annually-averaged polar values have 
dropped below 100 DU. In addition, this midlatitude maximum slowly shifts equatorward in both 
hemispheres. A substantial portion of the remaining total column is found in the troposphere (ap- 
proximately 20 DU in the Antarctic, 24 DU in the Arctic, and 20 DU in the tropics, with a globally- 
245 averaged tropospheric column of 27 DU). 

In the WORLD AVOIDED simulations polar total ozone also shows severe losses. Over the Arctic 
(Fig. 4a), ozone decreases from values near 500 DU in the 1960-1980 period to less than 100 DU 
by 2065 (an 80% depletion over approximately an 80-year time span). Approximately 24 DU of 
this remaining Arctic column is in the troposphere. The inset images show the April monthly aver- 
250 ages for 1980, 2020, and 2040. In 1980, the large column amounts cover the Arctic region (values 
> 500 DU). By 2020, a distinct “ozone hole” minimum has developed over the Arctic with a low 
value near the pole that is less than 200 DU. By 2060 values of less than 100 DU are found over 
the Arctic. The positive gradient of ozone between the midlatitudes and the pole seen in 1980 has 
disappeared by 2060. 

255 Over Antarctica (Fig. 4b), the October average drops below 100 DU in approximately 2025 and 
continues to slowly decrease to about 50 DU by 2060. Of this 50 DU, approximately 20 DU is found 
in the troposphere. The inset images show the October averages for 1980, 2020, and 2060. In 1980, 
a distinct ozone low is seen over Antarctica. By 2040, this low has considerably deepened, and by 
2060 extremely low values are observed across Antarctica and into midlatitudes. These October 
260 images also show the very large depletions of ozone in the midlatitude collar region. 

Antarctic ozone, at the profile peak (« 50 hPa), reaches 100% loss by the year 2000 (saturation), 
as observed in ozonesondes (Solomon et al., 2005). This saturation occurs because of the colder 
Antarctic temperatures and thereby greater coverage of PSCs and cold sulfate aerosols. The het- 
erogeneous reactions on the surfaces of these PSCs and aerosols fully activate chlorine, leading to 
265 massive ozone loss. Hence, in Fig. 4b we see a rapid ozone decline before 2000 with a slower de- 
cline after 2000 as ozone losses expand into the regions above and below 50hPa. In contrast, Arctic 
total ozone decreases linearly because saturation effects do not occur (e.g.,Tripathi et al., 2007). 

In addition to radical changes in the spring polar column amount in the WORLD AVOIDED sim- 
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ulation, the annual cycle of total ozone over the polar regions is also radically altered. Figure 5 
270 displays the column values versus the day of the year over the Arctic (Fig. 5a) and the Antarctic 
(Fig. 5b) at 10-year increments from 1980 to 2065. The Arctic (Fig. 5a) peak levels in the early 
years (1980-1990) occur in the springtime. Large depletions in the Arctic spring become apparent 
by the late- 1990s, and gradually worsen in later years. As seen in Fig. 4a, the April ozone reaches 
a 220-DU level in about 2030. The observed autumn through winter increase of ozone is clearly 
2 75 seen in 1980 (Bowman and Krueger, 1985), but by 2020 this has been inverted to a decrease. A 
comparison of the 1980 values with 2040 reveals a complete phase shift of the annual cycle. 

The Antarctic annual cycle (Fig. 5b) shows a worsening spring situation. The development of 
the ozone hole is clearly seen in the spring by comparison of 1980 to the 1990 and 2000 curves. 
The ozone hole continues to worsen into the 21st century with lower values and earlier onsets of 
280 the minimum value. Also note that the spring-to-summer increase of ozone has weakened by 2030, 
thereby creating a year-round ozone hole. This weak summer increase results from the decreased 
advection into the Antarctic region as ozone is depleted in both the upper stratosphere and midlati- 
tude middle and lower stratosphere. By 2050 the annual cycle of ozone over Antarctica is relatively 
small because of the very strong in situ depletion and the lack of advective resupply of ozone to the 
285 Antarctic region. 

The WORLD AVOIDED simulation has its largest losses in two distinct vertical layers: the lower 
stratosphere (200-30 hPa), and in the upper stratosphere (2 hPa). Since most of the ozone is found 
in the lower stratosphere, the lower stratospheric losses dominate the column losses seen in Fig. 2 
through Fig. 5. Figure 6 displays the 2060 difference in annually-averaged ozone losses as a function 
290 of altitude and latitude between the WORLD AVOIDED and the reference future simulations. The 
annually-averaged losses in the 200-30 hPa region of Antarctica exceed 90% in the region extending 
out to the edge of the polar vortex (the peak of the zonal mean winds is colocated with the polar 
vortex edge). Arctic lower stratospheric annually averaged losses also exceed 90% by 2060. In 
addition to the polar loss, a large tropical lower stratospheric loss (> 70%) is also observed in the 
295 70-30 hPa layer. Morgenstern et al. (2008) showed a very similar pattern of upper stratospheric 
losses and polar losses in comparison to our WORLD AVOIDED 2020 to 2000 differences. 

Upper stratospheric losses exceed 60% over a large latitude width in the 5-0.5 hPa region, with 
peak losses above 70% in the midlatitudes. Middle stratosphere losses are somewhat less (40- 
60%) in the 5-3 hPa layer. In addition to the stratospheric losses, ozone losses extend into the 
300 troposphere (the tropopause is indicated in magenta on Fig. 6). The GEOSCCM model relaxes to a 
fixed ozone climatology in the troposphere with a 5-day time scale, hence, the decrease of ozone in 
the troposphere is caused by advection of ozone-depleted air i nto the troposphere. The small increase 
of ozone just below the tropical tropopause is probably related to the increased vertical lifting. 

The ozone loss in the tropical lower stratosphere occurs extremely rapidly in the WORLD AVOIDED 
305 simulation in the six-year period from 2052 through 2058. Figure 7 shows the annual average of 
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ozone in the tropics (10°S-10°N) at 50 hPa from the WORLD AVOIDED simulation. Ozone shows 
a slow decline from the 1960s to the 2050 period as the residual circulation accelerates (Oman 
et al., 2008). The inset to Fig. 7 shows the February ozone values plotted against the corresponding 
temperature. The ozone decrease quickly accelerates as temperatures drop below about 200 K in 
310 the 2052-2058 period. This is a result of ozone loss cycles that are accelerated by heterogeneous 
processes in this period (to be discussed in the next section). 

5 Chemistry 

Stratospheric ozone levels in the WORLD AVOIDED simulation are severely depleted by the chlorine 
catalytic loss cycles. Furthermore, the relative fractions of chlorine, bromine, nitrogen, and hydrogen 
31 5 compounds are also drastically changed. The CIO + CIO and CIO ~F BrO catalytic cycles dominate 
ozone loss in the Arctic and Antarctic (Molina and Molina, 1974), while in the upper stratosphere 
the loss is dominated by the Cl + 0 3 and Br + 0 3 catalytic loss cycles. Figure 8 displays the Cl y 
partitioning versus year in the Arctic (80°-90°N average, 52hPa— the model level nearest 50hPa) 
for 1 March of each year. The Cl y does not grow smoothly (as shown in Fig. 1) because of year- 
320 to- year transport variability and because the Cl y is only shown for 1 March. The individual species 
are indicated by the color bar (dominant species are HC1, CIONO2, and CI2O2). Up to about 2025, 
the reactive chlorine provides a major fraction of the Cl y , but by 2035 the reactive species are only 
a minor component with HC1 dominating Cl y . This region is still largely in polar night (solar zenith 
angle > 87.7°), hence, the reactive Cl y is mainly found in CI2O2. The ozone levels are indicated 
325 in Fig. 8 by the superimposed black line. Polar lower stratospheric ozone falls below 1 ppmv in the 
2020s and remains low through 2065 with some minor variations that are the result of transport and 
temperature variations (less loss in warmer winters and more loss in colder winters). The low levels 
of ozone in these later years are accompanied by the preponderance of HC1 in the Cl y reservoir. This 
HC1 dominance results from the reaction of the Cl with CH4 and the disappearance of C10 x species 
330 when ozone levels approach zero (Douglass et al., 1995). 

A notable result of the WORLD AVOIDED simulation is the extremely rapid tropical lower strato- 
spheric depletion in the 2052-2058 period (see Fig.7). The inset of Fig.7 shows the strong relation- 
ship of temperature to lower stratospheric ozone, suggesting that the rapid ozone loss results from 
a nonlinear interaction of the cooling of the tropical lower stratosphere and heterogeneous chem- 
335 istry. Figure 9 displays the tropical lower stratospheric Cl y versus year partitioning at 52hPa and 
10°S-10°N for 1 February of each year. HC1 (80-85%) and C10N0 2 (10-15%) dominate Cl y in 
the 1975-2050 period. During the 2052-2058 period, the levels of reactive chlorine (C1 2 0 2 and 
CIO) suddenly increase from a few percent to about 20% of the total, while C10N0 2 drops from 
20% to only 1-2% of the total. Concurrent with the rapid increase of the reactive chlorine, the 
340 ozone drops to near zero (black lines in Fig.7 and Fig.9). This rapid repartitioning of the chlorine 
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budget results from the increased reactivity of HC1 and CIONO 2 on sulfate aerosol particles as the 
lower stratosphere cools. Such increased reactivity has been observed over Antarctica (Kawa, et ah, 
1997). “Ozone hole” chemistry appears in the tropical lower stratosphere in about 2052, leading to 
complete lower stratospheric ozone loss by 2058. 

345 6 Dynamics and Transport 

The large ozone depletions in the WORLD AVOIDED simulation lead to remarkable changes in 
dynamics and transport, including large temperature decreases by the middle of the 21st century. 
Figure 10 displays the temperature change between the WORLD AVOIDED and reference future 
simulations during the 2055-2065 decade. Cooling is observed everywhere in this simulation ex- 
350 cept in the troposphere and at 36 km over Antarctica. The cooling at 1 hPa exceeds 20 K. Shortwave 
heating at 1 hPa is directly proportional to the ozone concentration. As ozone decreases the short- 
wave heating decreases and the upper stratospheric temperature decreases. 

The WORLD AVOIDED cooling in the tropical lower stratosphere results from the combined ef- 
fects of an increase of the vertical lifting and decreased shortwave heating in the tropical lower strato- 
355 sphere. The shortwave heating remains roughly constant up to 2045 because of the balancing effect 
of the UV penetration and ozone changes. As ozone decreases in the upper stratosphere greater UV 
flux penetrates to 50 hPa, but as 50-hPa ozone decreases the UV flux absorption decreases, resulting 
in only very small changes to the shortwave heating. As the increasing tropical lower stratosphere 
lifting continues to cool the lower stratosphere, tropical lower stratospheric temperatures cool below 
360 the threshold for forming stratospheric clouds, and heterogeneous chemical processes lead to more 
ozone loss. The larger ozone losses then lead to less shortwave heating which further cools the lower 
stratosphere, producing a positive feedback that accelerates the ozone loss. This positive feedback 
leads to the complete collapse of lower stratospheric ozone between 2052 and 2058 (Fig. 7). This 
same ozone collapse is also simulated in the coupled 2-D model WORLD / VOIDED . 

365 Changes of ozone and temperature result in zonal-mean zonal wind changes. Figure 1 1 displays 
the change during January between the reference future and WORLD AVOIDED simulations aver- 
aged over the 2055-2065 period. In the Southern Hemisphere (SH) the mid-summer easterlies in the 
reference future simulation are replaced with westerlies extending to I hPa in the WORLD AVOIDED 
simulation. At 1 0 hPa and 65°S the change is approximately 45 ms" 1 , while the WORLD AVOIDED 
370 mean winds are slightly higher than 40 m s” 1 . The complete reversal of the usual summer easterlies 
indicates that the SH stratosphere has slipped into a permanent winter across the seasons by 2065. 
The Northern Hemisphere (NH) winter polar vortex is slightly weaker in mid- winter in the WORLD 
AVOIDED , but this difference is not statistically significant. In addition to these polar vortex results, 
the upper side of the subtropical jet (30° N, 70 hPa) is about 5ms” 1 stronger in January. 

The large ozone losses also induce changes in the transport of trace species in the stratosphere. 
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Figure 12 shows the difference in mean age-of-air between the WORLD AVOIDED and reference fu- 
ture simulations during the 2055-2065 period. The mean age is calculated in the model by advecting 
a tracer that increases linearly with time, and then computing the time difference at various points 
of the atmosphere with this age tracer at the tropical tropopause (100 hPa, 20°S~-20°N). Oman et al. 
380 (2008) have recently shown that there is near linear decrease of mean age over the late 20th and 21 st 

centuries in the reference past and reference future simulations. In the WORLD AVOIDED simu- 
lation, the large ozone depletion leads to additional decreases of the mean age-of-air by more than 
one year in the extratropical stratosphere. Much of this age decrease arises from the increase of the 
vertical lifting. The residual circulation is shown in Fig. 10 and Fig. 12 as the transparent white 
385 arrowed lines. Vertical lifting in the tropics (50hPa, 10°S-10°N) has changed from 15md _1 in the 
reference future to 25 m d“ 1 in the WORLD AVOIDED simulation (a 65% increase in upweiling). 
In addition, there is decreased rising motion above 50 hPa (point A in Fig. 12). This fresh young air 
from the troposphere then spills poleward decreasing the mean age in the extratropics (point B). 

Circulation changes are also apparent in the middle- to-upper stratosphere. The downward circu- 
390 lation (point C in Fig. 12) is caused by the changes in the mean zonal winds. As noted by (Stolarski 
et al., 2006), the ozone losses caused by the Antarctic ozone hole help the polar vortex persist into 
the summer period. The presence of westerly winds in the late spring and summer allows planetary 
waves to propagate vertically and deposit easterly momentum in the upper stratosphere, inducing a 
poleward and down ward circulation that would not normally be present. This increased poleward and 
395 downward circulation is exacerbated in the WORLD AVOIDED simulation by the additional ozone 
losses. The opposite problem is found in the NH (point D). The planetary wave deposition of easterly 
momentum is reduced in the upper stratosphere, reducing the poleward and downward circulation. 
The small horizontal and vertical gradients of mean age-of-air in the mid- to- upper stratosphere lead 
to only small changes of age in these regions. 

400 7 UV changes 

As a result of the large ozone depletions, the surface UV reaches extreme levels. To illustrate these 
changes we have calculated the surface UV flux in the northern midlatitudes at the height of sum- 
mer (2 July at local noon) for the WORLD AVOIDED simulation. UV calculations were done for 
a cloud-free atmosphere using Atmospheric Laboratory for Applications and Science 3 (ATLAS-3) 
405 extraterrestrial (ET) solar flux and the Total Ozone Mapping Spectrometer radiative transfer (RT) 
code (TOMRAD). In the TOMRAD code, the atmosphere is assumed to be plane-stratified, with 
Rayleigh scattering (Bates, 1984), WORLD AVOIDED ozone and temperature profiles, and pre- 
scribed surface albedo. The solution method is based on the successive iteration of the auxiliary 
equation of the RT (Dave, 1964, 1965; Dave and Furukawa, 1966). The calculations of atmospheric 
41 0 transmission were done at the original sampling wavelengths of the ozone cross-sections measured 
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by Bass and Paur (1985); Paur et al. (1985), with steps ^0.05nm. The transmittance values were 
linearly interpolated to the vacuum wavelengths of the high resolution ET solar flux data (ATLAS-3 
Solar Ultraviolet Spectral Irradiance Monitor ET) and multiplied by the ET spectrum. The TOM- 
RAD RT model compared well with spectral irradiance measurements and the output from other RT 
41 5 models run during the International Photolysis Frequency Measurement and Model Intercomparison 
campaign (Bais et al., 2003). 

Figure 13 displays the flux as a function of wavelength for 1980, 2040, and 2065. The vast 
majority of the UV increase occurs at wavelengths less than 310nm. The inset to Fig. 13 shows the 
ratio of the 2020, 2040, and 2065 UV flux to the 1980 UV flux. At wavelengths below 308 nm the 
420 UV flux has more than doubled by 2065 , and the flux has increased by a factor of 1 0 4 for wavelengths 
less than 288 nm. At 280 nm in 2065, the WORLD AVOIDED UV flux is 10 8 9 times stronger than 
the 1980 flux level. 

We have calculated the UV index for midlatitude NH conditions in mid-summer (a relatively 
densely populated zone on the globe with a number of large cities (e.g., Washington, D. C., U. S. A.). 
425 Figure 14 shows this UV index for the WORLD AVOIDED , reference future , and fixed chlorine 
simulations. For mid-summer clear-sky conditions, the UV index is normally very high at local 
noon; with a typical time to produce a perceptible sun burn (type II skin) of about 10-20 minutes. 
The mid-summer period has the highest UV index over the course of the year. The differences 
between the reference future and fixed chlorine simulations show increases on the order of 5-10% 
430 by the year 2000. The WORLD AVOIDED simulation has the UV index nearing 15 by 2040 and 
exceeding a value of 30 by 2065. This extreme value of the UV index would reduce the perceptible 
sunburn time from approximately 15 minutes to about 5 minutes. The largest UV index values are 
observed in the tropics where ozone columns are small, and the sun is directly overhead. 

These extreme UV increases would also lead to large increases of skin cancer (Slaper et ah, 1996). 
435 In the same manner as estimating the UV index, we apply a DNA damage action spectrum to our IJV 
spectrum calculations (Setlow, 1974). The DNA damaging UV increases by approximately 650% 
between 1980 and 2065. 

8 Summary 

In this study we have simulated the effects on the stratosphere of a steady growth of ozone depleting 
440 substances and compared those results to the normal expectations of the evolution of the stratosphere 
in the 21st century. As was suggested by the original Molina and Rowland (1974) study on CFCs in 
the stratosphere, large concentrations of ODSs in the atmosphere would have virtually destroyed the 
majority of the ozone layer by 2065 (global annual average losses > 60%). 

Very large ozone losses are computed at all latitudes in the stratosphere of the WORLD AVOIDED 
445 simulation. The largest losses in both % and DU are in the polar latitudes. However, surprisingly 
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large losses also occur in tropical latitudes as a result of heterogeneous chemical processes that 
occur in the 2052-2058 period. This “polar chemistry” in the tropics begin to appear in the lower 
stratosphere as a consequence of cooling resulting from increased vertical lifting. 

Most of the ozone depletion occurs in the lower stratosphere (12-24 km) and the upper strato- 
450 sphere (^42 km). Losses inside the polar vortex in both hemispheres exceed 90% by 2065. At 
20 km, annual average tropical levels are reduced by 80%. Ozone losses in the upper stratosphere 
exceed 60%. Middle stratosphere ozone losses are smaller, but still greater than 40%. 

As in current observations, the largest losses occur in the spring polar regions. Over Antarctica, 
the ozone hole develops quickly in the 1980-2000 period, with somewhat slower trends after 2000 
455 as a result of the complete ozone destruction in the lower stratosphere. October column values 
drop from 400 DU (somewhat higher than historic observations) to values of about 60 DU by 2065 
(85% loss). By 2040, the SH stratosphere ozone values are so low that interannual variability of the 
dynamics and transport has virtually no impact on total ozone, such that the total ozone internannual 
variability is near zero. In the Arctic, ozone declines are more linear with April total ozone dropping 
460 from values near 500 DU to about 100 DU in 2065 (an 80% decline). The annual ozone mixing 
ratios in the Antarctica and Arctic lower stratosphere (15-20 km) are nearly zero. 

The severe ozone depletions lead to interesting dynamical changes. The SH westerly circulation 
persists virtually year round by 2035, although temperatures still exceed 195 K in the late spring. 
This westerly circulation allows the vertical propagation of planetary scale waves (westerly winds 
465 are necessary for the vertical propagation of large scale Rossby waves) (Charney and Drazin, 1961). 
The deposition of easterly momentum from these waves then drives a poleward and downward cir- 
culation. In the NH stratosphere, the mid-winter dynamical changes are surprisingly small. The 
stratospheric polar night jet remains virtually the same. By 2055, the vertical circulation in the trop- 
ics (50hPa, 10°S-10°N) has increased by about 65%. This increased upwelling cools the lower 
470 stratosphere leading to the onset of heterogeneous chemistry and the collapse of ozone in the lower 
stratosphere in the 2052-2058 period. 

Our WORLD AVOIDED simulation is limited in a few respects. First, the simulation was per- 
formed with specified SSTs that were simulated by the NCAR COM3 using the IPCC Alb scenario. 
Hence, our tropospheric temperatures do not correctly respond to the increased CFG radiative fore- 
475 ing or stratospheric changes. Second, the model does not include tropospheric chemistry (i.e., the 
tropospheric ozone values are relaxed to the Logan (1999) climatology) and thus tropospheric ozone 
changes are not discussed in this study. Because tropospheric ozone is not explicitly modeled, the 
radiative forcing on the troposphere of tropospheric ozone is not correct. Since the tropospheric 
changes are not correctly simulated, the complete tropospheric effect on the stratosphere is not cor- 
480 rectly simulated. Third, since the tropospheric chemistry is not modeled, the simulated surface UV 
changes are limited to those only caused by the stratosphere. Fourth, the GEOSCCM is forced 
by scenario estimates of surface mixing ratios of both greenhouse gases and CFCs. The increased 
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stratospheric circulation will shorten the lifetimes of CFCs, The increased UV that penetrates to the 
troposphere might increase OH, increasing the oxidizing capacity of the troposphere and thereby 
485 shortening the lifetimes of gases such as CH 4 and HCFCs. Hence, feedback of the UV and circula- 
tion would affect the evolution of both climate and ozone depleting gases. 
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Table 1. GEOSCCM simulations description 


Simulation 

Year range 

ODS scenario 

Fixed SSTs 

Reference past 

1950-2004 

Ab a 

Observations: HadISSTl b 

Reference future 

1996-2099 

2000-2099 

Ab a 

A: HadGEMl c 

B: NCAR CCSM3 SRESA1B 

Fixed chlorine 

1960-2100 

Ab a , fixed to 1960 

1960-2000: Observations: HadISSTl b 
2001-2100: NCAR CCSM3 PCMDI 

WORLD AVOIDED 

1974-2065 

+3% per year 

1974-2049: NCAR CCSM2 SRESA1B 
1974-2049: NCAR CCSM3 SRESA1B 


a Montzka et al. (2003) 
b Rayner et al . (2003) 
c Johns, et al. (2006) 
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Fig. 1. EESC versus year for the world avoided as globally-averaged from the model’s inorganic chlorine 
and bromine at 4.5 hPa (thick black line). The magenta lines show scenarios from the Montreal Protocol, the 
London Amendments, and the Copenhagen Amendments. The red line indicates the estimate of ODS evolution 
in the Ab scenario and the blue line indicates observationally based EESC (Montzka et al., 2003). The green 
line shows the fixed level of 1960 chlorine (1.2ppbv). EEC1 for the WORLD AVOIDED is shown as the grey 
line. The time lag between the model EESC and the EEC1 is accounted for by the transit time from the Earth’s 
surface to the upper stratosphere (£3 4-6 years). 
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Fig. 2. Annually-averaged global ozone for the WORLD AVOIDED (solid black), reference future (r cd), fixed 
chlorine (green), and reference past (blue) simulations. The curves are smoothed with a Gaussian filter with a 
half-amplitude response of 20 years, except for the WORLD AVOIDED , which is unsmoothed. The dashed line 
shows the 2-D coupled model simulation of the “world avoided”. The grey-shaded inset shows the WORLD 
AVOIDED total ozone plotted against global annually-averaged EESC at 4.5 hPa from Fig. 1 . 
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Latitude 


Fig. 3. Annually -averaged total ozone versus latitude for various years from the WORLD AVOIDED (black), 
reference future (red ), fixed chlorine (green), and reference past (blue) simulations. The pole-to-pole area- 
weighted averages are shown in Fig. 2. 
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Fig. 4. Arctic total ozone in (a) April and (b) October for the WORLD AVOIDED (black), reference future 
(red), fixed chlorine (green), and reference past (blue) simulations. The curves are smoothed with a Gaussian 
filter with a half-amplitude response of 20 years, except for the WORLD AVOIDED , which is unsmoothed. The 
inset false-color images show (a) April and (b) October averages for 1980, 2020, and 2060 with 20-DU color 
increments (see inset scale). 
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Fig. 6. Percentage total ozone difference between the WORLD AVOIDED and the reference future simulations 
for a decadal average (2055 to 2065). The dashed contours are in 20% increments while colors are in 10% 
increments. The magenta line shows the WORLD AVOIDED tropopause. The white lines show the zonal-mean 
zonal winds (easterlies are dashed) for the WORLD AVOIDED simulation. 


26 


0 3 (ppmv) 



Fig. 7. Annually-averaged ozone (50hPa, 10°S-10°N) for the WORLD AVOIDED (black), reference future 
(red), fixed chlorine (green), and reference past (blue) simulations. The curves are smoothed with a Gaussian 
filter with a half-amplitude response of 20 years, except for the WORLD AVOIDED , which is unsmoothed. The 
inset figure shows the February WORLD AVOIDED ozone versus temperature (also at 50hPa, 10°S-10°N). 
The years in 5-year increments are highlighted in magenta. 
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Fig. 8. Arctic lower stratosphere sand chart of Cl y versus year from the WORLD AVOIDED simulation for 1 
March, 52hPa, and 80°~90°N. The individual contributions of various species to Cl y are shown in color (see 
color bar). Also shown is ozone averaged over the same region (black line). 
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Fig. 9. As in Fig. 8, except for 1 February and 10°S-10°N. 
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Fig. 10. Annually-averaged temperature difference (K) between the WORLD AVOIDED and the reference future 
simulations for the 2055-2065 period. Contour increments are 2 K, while color increments are 1 K. The thick 
black (red) line shows the WORLD AVOIDED ( reference future ) tropopause. The white lines show streamlines 
for the residual circulation differences between the WORLD AVOIDED and the reference future simulations 
over the period. 
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Fig. 12. As in Fig. 10 except for mean age-of-air differences in units of years. Contour increments are 0.2 
years, while color increments are 0.1 years. 
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Fig. 13. UV flux (W m“ 2 nm -1 ) is shown for the WORLD AVOIDED simulation in the years 1980, 2040, and 
2065 (thick black lines). The UV flux is calculated using the July 30°~50 C N zonal-mean ozone and temperature 
profiles, and assuming a time of local noon on 2 July. The thick grey line shows the erythemal (sunburn) action 
spectrum (right hand axis). The thin darker grey line shows UV flux at the top of the atmosphere. The inset 
figure shows the ratio of the UV flux for the years 2000, 2020, 2040, and 2065 to its value in 1980. The 
horizontal grey line in this inset indicates a doubling of the flux from 1980 conditions. 
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